Cold atom simulation of interacting relativistic quantum field theories 
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We demonstrate that Dirac fermions self-interacting or coupled to dynamic scalar fields can emerge 
in the low energy sector of designed bosonic and fermionic cold atom systems. We illustrate this 
with two examples defined in two spacetime dimensions. The first one is the self-interacting Thirring 
model. The second one is a model of Dirac fermions coupled to a dynamic scalar field that gives 
rise to the Gross-Neveu model. The proposed cold atom experiments can be used to probe spectral 
or correlation properties of interacting quantum field theories thereby presenting an alternative to 
lattice gauge theory simulations. 
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Introduction:- To date there exist several important 
open problems in quantum field theories (QFTs) rang- 
ing from the convergence of the scattering matrix to the 
infrared behaviour of quantum chromodynamics (QCD). 
Such problems can be probed analytically only by non- 
perturbative methods that seem to be intractable in three 
and four dimensions. As an alternative, numerical tech- 
niques or quantum simulations can play a central role in 
obtaining insight into the Standard Model. 

Recently, much interest is focused on simulating QFTs 
with cold atoms [TH5]. In the relativistic domain, 
these studies are usually limited to generating Dirac 
fermions [H [7] and background fields [SHIO] . In this let- 
ter, we propose the realisation of Dirac fermions inter- 
acting with dynamic fields by employing cold atoms in 
optical lattices. This constitutes the first step towards 
the quantum simulation of a general model with cou- 
pled relativistic fermionic and bosonic fields. Cold atoms 
provide a controlled medium with well understood in- 
teractions. A wide range of quantum optics and atomic 
physics technology allows for the preparation, manipula- 
tion and detection of a variety of interesting many body 
phenomena. Employing cold atoms gives us the possibil- 
ity to consider two, three and four spacetime dimensions, 
to tune the couplings of the interactions and to explore 
the behavior of multi-component fields. Compared to nu- 
merical simulations of lattice gauge theories on comput- 
ers, a physical simulation on a quantum system naturally 
overcomes the sign problem 

Here we show how to simulate a two-dimensional 
self-interacting model of Dirac fermions, known as the 
Thirring model [12l [13] and two-dimensional Dirac 
fermions coupled to a scalar field that is equivalent to 
the Gross-Neveu model. The Hamiltonians of these sys- 
tems are supported on one spatial dimension. The neces- 
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sary building blocks are the Dirac Hamiltonian, which 
describes relativistic fermions, and the interaction of 
fermions with themselves or with a dynamic scalar field. 
This goes beyond previous proposals concerned with dy- 
namical fermions coupled to classical fields. We show how 
the required components naturally emerge in the low en- 
ergy sector of specifically designed lattice Hamiltonians. 
The Dirac operator describes the continuum limit of cer- 
tain fermionic lattices [SllMllTS], as in graphene |T6], that 
remarkably give rise to both the spin and the linear dis- 
persion relation. Compared to previous approaches [S] 
here we realise the Dirac operator only by the free tun- 
nelling of single species of atoms. Slow spatial varia- 
tions of the lattice couplings result, in the Dirac pic- 
ture, to an interaction with a scalar field background jlO) . 
When these distortions are caused by coupling the lattice 
fermions to bosons, the resulting fields become dynamic. 
Self-energy terms can be implemented giving rise to a 
variety of interesting QFTs. The presented models are 
exactly solvable and serve for demonstrating the ability 
to simulate important properties of the Standard Model 
such as dynamical symmetry breaking and mass genera- 
tion with cold atoms. 

Thirring model:- Our starting point is the two- 
dimensional Thirring model |121 113] . It describes inter- 
acting fermions with the Hamiltonian 



= dx 
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Here the 7's are two-dimensional matrices satisfying 
= "^Viiv for /X, 1/ = 0, 1 with 77^^ = diag(l,-l), 
75 = 7o7i E^nd ^I' is a two-dimensional spinor with 
^ = ^-t^Q. The mass of the fermions is mg, g is their 
dimensionless self-interaction coupling constant and Vg is 
the sound velocity taken in high energy to be the speed 
of light. This model has exciting physics with the mass- 
less case, mo = 0, being equivalent to free bosons and 
the massive case being equivalent to the Sine-Gordon 
model [T7]. 




FIG. 1: (a) The one-dimensional superlattice with tunnelhng 
fermionic atoms that simulates Dirac fermions. Each unit cell 
includes two fermion sites, a and b. An alternating distortion 
of the tunnelling couplings J and J — 5 gives rise to the mass 
term, S. (b) The energy dispersion relation as a function of 
momentum, p. At half filling and for 5 — the low energy be- 
haviour is linear with respect to p, allowing the Dirac operator 
description. 



To simulate the Thirring model in an atomic system 
we consider a linear bicolourable fermionic lattice with 
spacing I (see Fig. [TJa)) which is subject to the Hubbard 
Hamiltonian 



H 



27r ^ [ - J{alb, + bla,+i) + Sajbi + H.c. 



+Ua]a^blb, . (2) 



Here J is the tunnelling coupling between neighbouring 
sites of the lattice of the same fermionic atoms a and b 
with {aijOj} = {bi,b^j} = Sij and all the other anticom- 
mutators vanishing. The a, b index is a spatial distinction 
within the unit cell that allows the encoding of the spin 
degree of freedom. The tunnelling distortion S occurs on 
alternating links, as shown in Fig. [ija), while U is the 
interactions coupling between fermions in the same cell. 

If we diagonalise the J-term of the Hamiltonian in 
the a, b basis we find the dispersion relation E± (p) — 
±2 1 cos y I plotted in Fig. [l|b). It can be easily seen that 
there is a single Fermi point, P = w/l, for which E± (P) = 
0. If the lattice is half filled with fermions, which occupies 
the valence band completely, the behaviour of the small 
energy fluctuations is governed by the Hamiltonian ex- 
panded around P. Setting p = P -I- A: for <C 1// we 
obtain a dispersion that is linear in momentum, k. Hence, 
the fermionic tunnelling term around the Fermi point can 
be efficiently described by the relativistic Dirac Hamil- 
tonian Jl J dkipl.a2k'ipk, where ip — (a? corresponds 
to the one-dimensional version of the Kogut-Susskind 
fermions |14l I15j . To assign the appropriate dimensions 
to the Dirac fermions we set ip — \/2l'i'. In this way we 
obtain the massless free Dirac fermion term of ([T]) written 
in the momentum representation with 70 = cri, 71 = ia^ 



FIG. 2: The regularised mass A/, in units of h/{vj), as a 
function of the tunnelling disorder S for various interaction 
strengths U. When no interactions are present the mass in- 
creases, as expected, linearly as a function of 5. The presence 
of interactions dramatically changes this behaviour even for 
moderate ratios of U /J. 



and = 2Tr ■ 21J. It can be verified that the continuum 
limit of ^ gives also rise to the rest of the terms in ([T]) 
with mov^/h = 27r • 2(5 and g = 2tt ■ 2Ul. The continuum 
limit corresponds to small lattice spacing. This is equiv- 
alent to restricting to the low energy sector of the system 
where states have a large wavelength support. These are 
exactly the states we are interested in for probing the 
infrared behaviour of QFT, such as the ground state and 
its gapped or gapless nature. 

The Hamiltonian ([2| can be realised with cold atoms as 
follows. The one-dimensional fermionic tunnelling term 
of Q appears when a fermionic gas is placed in an optical 
lattice with very tight confinement in the other two direc- 
tions. The (5-term corresponds to a uniformly decreased 
tunnelling coupling between sites of the same cell. It can 
be generated, e.g. by employing superlattices as seen in 
Fig. [ija) . This alternatively signifies that inhomogene- 
ity in the tunnelling coupling due to experimental im- 
perfections will generate a mass term as observed in [T]. 
The final U-teim results from the interaction between the 
atoms a and 6 present in the same cell [See supplemen- 
tary material at [URL will be inserted by AIP]]. 

The Thirring model is the simplest relativistic inter- 
acting QFT that one could implement in the laboratory 
with present technology. Its realisation can demonstrate 
the renormalisation of mass due to interactions. Indeed, 
mo 7^ is the fermionic mass of the classical theory. 
When the interactions g are introduced then a regulari- 
sation condition needs to be adopted, \p\ < A, where A 
is a momentum cutoff that excludes modes with unphys- 
ically high energy. Our system is naturally regularised 
due to the underlining lattice structure, where the mo- 
mentum cutoff is related to the optical lattice spacing by, 
A = tt/I. a quantum field theory is called renormalisable 
when this cutoff can be absorbed in the initial, bare pa- 
rameters of the model, such as the mass. The massive 
Thirring model is such a renormalisable theory where the 
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interactions g give rise to the following regularised mass 



S tan " Yrr U 1 TT J 

The corresponding coupling regime < U/ J < n is well 
within the realm of current experiments |18) . This ef- 
fective mass is exactly the energy gap above the ground 
state obtained for zero momentum excitations, while non- 
zero momenta give a continuum spectrum above the gap. 
Fixing the lattice spacing, I, according to the wavelength 
of the optical lattice, one can plot the effective mass as 
a function of the atomic parameters S and f/, as seen 
in Figure [2j Current experiments routinely probe such 
excitation gaps in cold atom systems [T]. 

The atomic simulation relies on direct exchange of 
atoms, while the value of M increases when the cou- 
plings U or 6 increase, in a non-perturbative way, or 
when the lattice spacing I decreases. The size of the gap 
can, hence, be efficiently controlled with a variety of ex- 
perimental means which alleviates the low temperature 
requirements jl9j . This facilitates the experimental mea- 
surement of the gap by a spectral analysis of the atomic 
system [1] . The renormalisation step would then necessi- 
tate to take the 'bare' mass mg decrease as a function of 
A so that a finite constant value of the 'physical' mass is 
obtained. It is intriguing that this renormalisation pro- 
cedure can be established experimentally by studying the 
spectral behaviour of the fcrmionic lattice system. Ob- 
serving such a strong renormalisation of parameters (see 
Fig. [2]) in a system of cold atomic gases provides a unique 
fingerprint of strong correlations. 

Fermion-scalar interaction and the Gross- 
Neveu model:- Next we consider a two-dimensional 
model where an A^-colour massless Dirac fermion 
n = 1,...,N interacts with a massive quantised scalar 
field $ according to the Hamiltonian 



IT 



dx 



(4) 



Here we assume summation over the colour index n, g is 
the coupling strength between the bosonic and fermionic 
fields and m is a mass scale that can be absorbed in $. It 
can be shown that this Hamiltonian yields exactly equiv- 
alent fermionic behavior as the Gross-Neveu model [2U] 
given by 



GN 



dx 



(5) 



The Hamiltonian i?GN describes massless fermions sub- 
ject to attractive interactions. This attraction causes the 
fermions to create bosonic pairs. Eventually, these com- 
posite bosons condense breaking spontaneously the Z2 
symmetry, 75 of Hamiltonian ([S]), thereby caus- 

ing the fermions to dynamically acquire mass. This fas- 
cinating property is very similar to the behaviour of the 



BCS theory of superconductivity or of four-dimensional 
QCD. Note that Q does not contain a kinetic term for 
the scalar field. This corresponds to a Yukawa theory 
with interactions mediated by infinitely massive fields, 
making their propagation point-like and resulting in the 
Gross-Neveu effective four-fermion interaction. 

We will now consider a cold atom system that gives 
rise to [2Tj, which would make it possible to observe 
the dynamical mass generation experimentally. The ki- 
netic term of the Dirac fermions, 5'„, can be produced 
by the same fermionic tunnelling term as in ([2]). In gen- 
eral, a variety of interaction terms can be generated be- 
tween bosonic and fermionic atoms. To conform with 
Q, we specifically want the bosonic modes, to cou- 
ple linearly to the fermionic ones, as dictated by the 
minimal coupling prescription. Such an interaction can 
result similarly to the mo-term of the Thirring model. A 
position dependent tunnelling distortion gives rise to a 
classical scalar field configuration. Formally, the quan- 
tisation of the scalar field is obtained by writing it as 
$ = (d^ + d) /\/2, where d is a bosonic mode. Substitut- 
ing it into gives the quantised Dirac fermion-scalar 
model. In the cold atom setting, this can be achieved 
by employing a bosonic condensate interacting with the 
lattice fermions as we shall see below. 

Consider a one-dimensional fermionic lattice super- 
posed with a one-dimensional bosonic lattice, as seen in 
Fig. [Sj We assume that the dynamics of the atoms is 
described by 

y = 27r ^ I [ - {J'al^^bna + Jbl iana+i) - 

i 

5a\aia^^ J)n^i ^'R.c] +Ua\ -f- (6) 

Here the a„ 's and 5„ 's are N different species of fermionic 
atoms (summation over n is assumed) and ai — D + 
di is an atomic condensate with particle density D and 
bosonic operator d. When the couplings of Q and the 
condensate density are appropriately tuned, then the low 
energy behaviour of the Hamiltonian H reproduces 
Indeed, for J' = J—SD^ and /i = —2D^U one reproduces 
the desired low energy behavior with g^ = 2tt-2S'^1/U and 
Vs = 2TT ■ 2Jl. To suppress spurious terms we take the 
weak fluctuation limit {d^d) <C D^. 

The terms inj6]) can be realised in the lattice config- 
uration of Fig. |3j by employing one-dimensional optical 
lattices filled with the appropriate species of atoms. All 
the required interactions naturally appear in cold atom 
settings and can be tuned, e.g. by Feshbach resonances. 
The J- and J'-terms of ([6| result from the tunnelling 
of the fermionic atoms along their lattice and contribute 
to the free Dirac fermion propagation. The i5-term re- 
sults from the interaction between bosonic and fermionic 
atoms that generate a fermionic tunnelling controlled by 
the bosonic populations. It gives rise to the interaction 
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Fermionic 
optical lattice 



■ Bosonic 
optical lattice 



FIG. 3: The one-dimensional optical lattice with tunnelling 
fermions and bosons that simulates the Dirac fermion-scalar 
field model. The bosonic sites are placed in-between the 
fermionic ones with double spacing. In this way the bosonic 
population on site i controls the fermionic tunnelling within 
the same cell. 

between the scalar and the fermion fields. In ^ there 
is no tunnelling term for the bosonic atoms. Hence, the 
required system comprises of independent BECs at each 
site of a bosonic lattice [22j[23]. The ?7-term describes the 
interaction between bosonic atoms on the same site and 
the /i-term is given by their chemical potential. Impor- 
tantly, small errors in the values of all of these couplings 
result in modifications of the parameters of the effective 
Hamiltonian ^ or the generation of similar terms [lO] . 
The persistence of the effective Dirac description in the 
low energy limit is a characteristic also present in the 
Thirring model simulation. 

It has been shown [20] tbat the initially massless 
fermions of the Gross-Neveu model dynamically acquire 
an effective mass. In terms of the atomic parameters, the 
emerging mass gap is given by 



TT TT JU 



(7) 



This result is exact in the large N limit but is expected 
to be a good approximation even for moderate, exper- 
imentally feasible, values of N of the order of 2 or 3. 
While the mass, M, goes to infinity when / goes to zero, 
the renormalisation prescription absorbs the infinity in 
the 'bare' coupling g. In the case of the optical lattice 
realisation we do not need to renormalise the coupling g 
as we are working with a fixed lattice spacing I. Such a 
simulation can detect the dynamical generation of mass 
and verify the predicted behaviour ([t]) as a function of 
the atomic couplings. 

Conclusions:- The above method can be straightfor- 
wardly generalised to other more complex QFTs. For 
example, to simulate quantum electrodynamics in two 
spacetime dimensions (QED2) we need a similar atomic 
setting to the fermion-scalar model described above. 
Apart from the addition of a kinetic term for the bosonic 
field we need to impose the Gauss constraint in the 
atomic level. An atomic model that can give rise to 
QED2 as well as its generalisation to four dimensions 
will be presented elsewhere. 

Going beyond QED one could realise the SU(A^) Yang- 
Mills gauge field coupled to an iV-colour Dirac fermion 
(see Q) that gives rise to QCD. The challenge faced 



at this point is to realise the higher order bosonic self- 
interaction terms that are necessary to simulate the non- 
linear behaviour of the Yang-Mills theory. It is a fasci- 
nating perspective to probe the generation of a mass gap 
in such a QCD simulation. In parallel, a wealth of possi- 
bilities opens for the simulation of supersymmetric QFTs 
or combining QFT and relativity in an atomic system. It 
is highly plausible that the experimental realisation of 
the simple models presented here provide a unique prob- 
ing tool into the open questions of interacting relativistic 
QFTs. 
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